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> 1 INTRODUCTION 



ABSTRACT 

We present mass models of a sample of 19 low surface brightness (LSB) galaxies and 
compare the properties of their constituent mass components with those of a sample 
of high surface brightness (HSB) galaxies. We find that LSB galaxies are dark matter 
dominated. Their halo parameters are only slightly affected by assumptions on stellar 
mass-to-light ratios. Comparing LSB and HSB galaxies we find that mass models 
derived using the maximum disk hypothesis result in the disks of LSB galaxies having 
systematically higher stellar mass-to-light ratios than HSB galaxies of similar rotation 
velocity. This is inconsistent with all other available evidence on the evolution of LSB 
galaxies. We argue therefore that the maximum disk hypothesis does not provide a 
representative description of the LSB galaxies and their evolution. Mass models with 
stellar mass-to-light ratios determined by the colors and stellar velocity dispersions of 
galactic disks imply that LSB galaxies have dark matter halos that are more extended 
and less dense than those of HSB galaxies. Surface brightness is thus related to the halo 
properties. LSB galaxies are slowly evolving, low density and dark matter dominated 
galaxies. 

Key words: dark matter — galaxies: kinematics and dynamics — galaxies: spiral 
— galaxies: fundamental parameters — galaxies: halos 



The discrepancy between the amount of matter implied by 
the Hi rotation curves of spiral galaxies and the amount of 
matter actually observed (in the form of stars and gas) is 
usually interpreted as dark matter (DM) halos that surround 
the directly observable parts of galaxies, although alter- 
native theories of modified Newtonian dynamics (Milgrom 
1983) also provide an efficient description of the observed ro- 
tation curves (Sanders 1996). Early studies by e.g. Kalnajs 
(1983) claimed that, based on the optical rotation curves 
that were used, there was no need to invoke DM within 
the optical radius for some galaxies. Hi observations, how- 
ever, showed a dramatically different picture (Bosma 1978, 
Begeman 1987) . Substantial amounts of DM were needed to 
describe the observed fiat rotation curves outside the optical 
disk. 

Measurements of the distribution of the DM must be ex- 
tracted from the rotation curve. This is usually done by com- 
puting the rotation curves of the visible matter and subtract- 
ing these from the observed total rotation curve. The resid- 
uals show the dynamical signature of DM. This procedure 
typically has the following free parameters: a length scale 
and density for the halo and a mass-to-light ratio (M/L)* 



of the stellar component. A major problem is that the value 
of {M/L)i, is not known a priori, and may differ from galaxy 
to galaxy. Many different values of {M/L)i, and therefore 
many combinations of halo parameters, yield equally good 
descriptions of the data (van Albada et al. 1985). It is for 
many galaxies possible to make a good fit to the rotation 
curve by completely ignoring the contribution of the stellar 
disk, even though it obviously plays an important role in 
the inner parts of many galaxies. Without additional knowl- 
edge to constrain (M/L)* (e.g. information about vertical 
stellar velocity dispersion or derived from stellar population 
synthesis models) it is not possible to unambiguously de- 
termine its value. Goodness-of-fit estimators are sensitive to 
small variations and uncertainties in the data (even when 
the data is of superior quality) (Lake & Feinswog 1989) and 
are thus less suited for simultaneously determining all free 
parameters in a fit. The mathematically "best" fit is there- 
fore not always the physically most meaningful one. 

Using the observation that the rotation curves derived 
for the stellar components of HSB galaxies can usually be 
scaled so that they almost completely describe the observed 
total rotation curves in the optical disks of these HSB galax- 
ies (e.g. Kent 1986), van Albada & Sancisi (1986) introduced 
the maximum-disk hypothesis. This hypothesis reduces the 
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number of free parameters in rotation curve fits by maxi- 
mizing the contribution of the stellar disk [i.e. (M/L)*] and 
thus minimizing the amount of DM invoked to explain the 
observed rotation curves. 

It is still unclear whether the maximum disk hypothesis 
is a realistic one. For some HSB galax;ies stellar population 
synthesis models yield values of (M/L)j, that are consistent 
with the maximum disk measurements (van Albada & San- 
cisi 1986). The presence of disk features such as spiral arms 
or bars argues for a high (M/L)^, (Athanassoula, Bosma 
& Papaioannou 1987). Kuijken & Gilmore (1989) showed, 
however, that there is only enough matter in the disk of 
our Galaxy for it to be about half-maximum disk. Mea- 
surements of the stellar velocity dispersions in HSB galaxies 
(Bottema 1995) show that maximum disk may overestimate 
the amount of matter present in the disk by some 60 per 
cent. 

As the maximum disk hypothesis gives for each galaxy 
the maramum possible (Af/L)*, any conclusions about struc- 
tural parameters of halos are only "lower limits." Neverthe- 
less, attempts have been made to find systematic relations 
between the observable properties of spiral galaxies and their 
dark halo properties. 

A general conclusion is that the importance of DM in- 
creases towards later types (Persic & Salluci 1991). Coupled 
with the finding that the surface brightness also decreases 
towards later types (de Jong 1995), this means that low 
surface brightness (LSB) galaxies arc a prime example of 
galax;ies that are extremely interesting for deriving struc- 
tural properties of DM halos. 

Those galaxies arc very late-type spirals with central 
surface brightnesses much lower than those of "classical" 
late-type HSB galaxies (de Blok, McGaugh & van dcr Hulst 
1996 [hereafter BMH96], Zwaan et al. 1995, van der Hulst 
et al. 1993; see also the reviews by Impy & Bothun 1997 
and Bothun, Impoy & McGaugh 1997). Their blue col- 
ors (de Blok, van dcr Hulst & Bothun 1995, McGaugh & 
Bothun 1994, Ronnback & Bergvall 1994), low star forma- 
tion rates (McGaugh 1992, van der Hulst et al. 1993, van den 
Hoek et al. 1997), low oxygen abundances (McGaugli 1994, 
Ronnback & Bergvall 1995) and high gas-fractions (BMH96, 
McGaugh & de Blok 1997), all give rise to the picture that 
these galaxies are slowly evolving and still in an early evo- 
lutionary state. 

Recently a number of rotation curves of LSB galax;ies 
have boon measured (BMH96). A preliminary investigation 
showed a systematic trend of increasing total mass-to-light 
ratio within a number of scale lengths with decreasing sur- 
face brightness, from values of ~ 1 for HSB galaxies to ~ 10 
for the dimmest measured LSB galaxies (Fig. 13 in BMH96), 
consistent with predictions made using the TuUy-Fisher re- 
lation (Zwaan et al. 1995). The increasing importance of 
DM makes an investigation into the halo parameters extra 
interesting as the parameters will depend loss strongly on 
assumed values of (M/L)*, and one can therefore be quite 
sure that the derived parameters are indeed close to the in- 
trinsic parameters of the halo. 

A comparison of a HSB and a LSB galaxy at identical 
positions on the TuUy-Fishor relation (i.e. at fixed V^max) (de 
Blok & McGaugh 1996) showed that LSB galaxy UGC 128 is 
almost a factor 10 less dense than the HSB galaxy NGC 2403 
(analogous with their low stellar and Hi surface densities). 



UGC 128 also has a factor of 7 higher total M / L (measured 
within a fixed number of scale lengths) and must therefore 
be more massive than HSB galaxy NGC 2403 (which has 
the same luminosity) . 

In this paper we will discuss disk-halo decompositions 
of the sample of LSB galaxies from BMH96 and compare the 
structural parameters of the dark and luminous mass com- 
ponents of HSB and LSB galaxies. In Section 2 we describe 
the sample; Section 3 covers various aspects of the effects 
on resolution and beam-smearing on the rotation curves. In 
Section 4 the actual decompositions are presented, while in 
Section 5 we discuss the results and compare them with what 
is known about other galaxies. In Section 6 the systematic 
relations between the masses of the difi^ercnt components 
(gas, stellar, dark) are discussed, while in Section 7 we dis- 
cuss the trends of volume and surface density as a function 
of surface brightness. In Section 8 our main results are sum- 
marized. In Appendix A a summary of the Bottema disk 
(Sect. 4.4) is given. 



2 THE LSB SAMPLE 

We refer to BMH96 for an extensive sample description. We 
have supplemented that sample with the rotation curves of 
an additional 5 LSB galaxies as presented in van der Hulst 
et al. 1993. The rotation curves in general have maximum 
rotation velocities between 50 and 120 km s^^, and rise 
more slowly than HSB curves with similar mgiximum ve- 
locity. Only a few of the rotation curves show clear signs of 
flattening towards the outermost radii. Optical properties 
of the sample are presented in McGaugh & Bothun (1994), 
de Blok et al. (1995), and van dcr Hulst et al. (1993). The 
properties of the sample galaxies are summarized in the top 
panel of Table 1. We adopt a Hubble constant of 75 km s~^ 
Mpc-^ 

As described in BMH96 the rotation curves were de- 
rived from major-axis position-velocity diagrams. Inclina- 
tions were determined from a comparison of optical, Hi and 
kinematical inclinations. The position angle and inclination 
wore both kept fixed during the derivation of the curves. 
Inspection of the position-velocity diagrams in Fig. 2 of 
BMH96 shows that not all rotation curves are well suited 
for a disk-halo analysis. A number of observations clearly 
suffer from low resolution or low inclination. 

In order to select only those galaxies where a disk-halo 
decomposition is a sensible exercise, we restricted the sample 
by rejecting those galaxies that did not meet one of the 
following criteria: 

(i) inclination i > 25°, 

(ii) radius R> 2 beams, 

(iii) asymmetries between approaching and receding sides 

of rotation curves < 10 km s~^. 

This resulted in the rejection of seven galaxies. These are 
listed in Table 2 with the reason for rejection. 

The choice for a 25 degree cut-off may seem a rather 
liberal choice, which in a sense we are forced to apply due 
to the selection biases against high-inclination LSB galaxies 
(Davies 1990, Schombert et al. 1992, Dalcanton & Shectman 
1996). An analysis of the scatter in the LSB TF relation 
presented in Zwaan et al. (1995) does show however that 
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Table 1. Properties of LSB and HSB samples 



Name 


D 




n 


ID — V 
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D 
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T/ 




(Mpc) 


(nr, o <T / 1 — I" ^ 

(mag/U ) 


(.Kpcj 


(mag) 


(mag) 


(kpcj 


(Kpc) 


(km s ^) 


LSB sample (this 


work) 














F561-1 


63 


23.3 


3.6 


0.55 


-17.8 


6.4 


10.1 


52 


F563-1 


45 


23.6 


2.8 


0.64 


-17.3 


5.0 


17.7 


111 


F563-V1 


51 


24.3 


2.4 


0.59 


-16.3 


3.2 


7.4 


30 


F563-V2 


61 


22.1 


2.1 


0.51 


-18.2 


6.8 


9.2 


111 


F564-V3 


6 


24.1 


0.4 


0.56 


-11.1 


0.3 


~2 


~40 


F565-V2 


48 


24.7 


2.7 


0.51 


-15.4 


2.3 


8.4 


51 


F567-2 


75 


24.4 


5.7 


0.61 


-17.4 


3.2 


11.3 


64 


F568-1 


85 


23.8 


5.3 


0.58 


-18.1 


7.4 


14.9 


119 


F568-3 


77 


23.1 


4.0 


0.61 


-18.3 


7.9 


16.5 


120 


F568-V1 


80 


23.3 


3.2 


0.57 


-17.9 


6.2 


19.0 


124 


F571-8 


48 


23.9"'= 


5.2= 


* 


-17.6" 


7.7 


15.6 


133 


F571-V1 


79 


24.0 


3.2 


0.55 


-17.0 


3.8 


14.6 


73 


F571-V2 


16 


* 


* 


* 


* 


* 


3.7 


45 


F574-1 


96 


23.3" 


4.3 


* 


-18.4" 


17.7 


15.4 


100 


F574-2 


88 


24.4 


6.0 


0.59 


-17.6 


3.8 


10.7 


40 


F577-V1 


80 


24.0 


4.3 


0.40 


-18.2 


6.6 


8.9 


30 


F579-V1 


85 


22.8" 


5.1 


* 


-18.8" 


11.1 


17.3 


100 


F583-1 


32 


24.1 


1.6 


* 


-16.5 


3.3 


14.6 


85 


F583-4 


49 


23.8" 


O '7 


* 


1 a no 
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52 
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nn 
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U5750 


56 


23.5" 


O.D 
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-10. / 


y.o 


ZL.O 


i 
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45 
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sample (Broeils 1992) 














DD0154 


4 


23.2 


0.5 


0.29 


-13.8 


1.0 


7.6 


48 


DD0168 


3.5 


23.4 


0.9 


0.22 


-15.2 


1.7 


3.4 


55 


DDO170 


12 


* 


1.3 


* 


-14.5 


2.5 


9.6 


66 


N55 


1.6 


21.3 


1.6 


0.38 


-18.6 


8.9 


10.2 


87 


N247 


2.5 


23.4 


2.9 


0.45 


-18.0 


7.3 


9.9 


108 


N300 


1.8 


22.2 


2.1 


0.56 


-17.8 


5.1 


10.6 


97 


N801 


79.2 


21.9 


12 


0.66 


-21.7 


36.4 


58.7 


222 


N1003 


11.8 


21.7 


1.9 


0.40 


-19.2 


9.2 


31.3 


115 


N1560 


3 


23.2 


1.3 


0.43 


-15.9 


4.3 


8.3 


79 


N2403 


3.3 


21.4'' 


2.1 


0.38 


-19.3 


8.4 


19.5 


136 


N2841 


18 


21.1'' 


4.6 


0.80 


-21.7 


22.6 


81.1 


323 


N2903 


6.4 


20.5'' 


2.0 


0.59 


-21.0 


11.7 


24.2 


201 


N2998 


67.4 


20.3 


5.4 


0.48 


-21.9 


28.2 


46.6 


214 


N3109 


1.7 


23.2 


1.6 




-16.8 


4.3 


8.2 


67 


N3198 


9.4 


21.6** 


2.6 


0.46 


-19.4 


11.4 


29.9 


157 


N5033 


11.9 


23.0'' 


5.8 


0.48 


-20.2 


18.3 


35.4 


222 


N5533 


55.8 


23.0 


11.4 


0.80 


-21.4 


26.2 


74.4 


273 


N5585 


6.2 


21.9 


1.4 


0.42 


-17.5 


4.7 


9.6 


92 


N6503 


5.9 


21.9 


1.7 


0.56 


-18.7 


5.4 


22.2 


121 


N6674 


49.3 


22.5 


8.3 


0.61 


-21.6 


29.8 


64.5 


266 


N7331 


14.9 


21.5** 


4.5 


0.70 


-21.4 


23.4 


36.7 


241 


U2259 


9.8 


22.3'' 


1.3 


* 


-17.0 


3.7 


7.6 


90 


U2885 


78.7 


22.0'' 


13 


* 


-22.8 


62.8 


72.5 


298 



a: converted from i?-band measurement 

b: derived from r-band profiles, converted using B = r + 1 (Kent 1986) 
c: /i()(face on) of F571-8 was computed assuming /i(facc on) = (zo//i)/i(edgc on) (i.e. 
optically thin), hiia.ce on) was found to be ~ 1.4/i(edge on) following method of 
Peletier et al. (1995). 
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Table 2. Rejected galaxies 



Name 



Limitation 



F561-1 R<2 beams; i < 25° 

F563-V1 R < 2 beams; gross asymmetry 

F564-V3 R<2 beams; low S/N 

F567-2 gross asymmetry; i < 25° 

F574-2 R < 2 beams; V(R) steeply rising 

F577-V1 R < 2 beams; gross asymmetry 

F579-V1 gross asymmetry 



this scatter does not increase significantly when increasingly 
lower-inclination galaxies are included until a value of 25 
degrees. 

The effects of lower resolution are more difficult to 
model, and will be described extensively in the next section. 



3 BEAM-SMEARING 

Early high-resolution studies of spiral galaxy kinematics us- 
ing synthesis radio telescopes did mostly target galaxies with 
large angular sizes (e.g. Bosma 1978, Begeman 1987, Broeils 
1992). The rotation curves of these galaxies are in gen- 
eral well-determined. Unfortunately there are only a limited 
number of galaxies in the sky that have suitably large an- 
gular sizes for such detailed investigations. Galaxies with 
smaller apparent sizes (because of smaller intrinsic sizes 
and/or larger distances) can only be observed at a smaller 
linear resolution. 

A side-effect of this lower resolution is that the line 
profiles will be artificially broadened and any sudden change 
in the velocity gradient within the beam will be smoothed 
out. The observed change in the velocity gradient will thus 
appear to be smaller and the rotation curve will appear to 
rise less steeply. This effect is called "beam smearing." 

Bosma (1978) investigated the behavior of a steeply ris- 
ing model rotation curve at different resolutions. He showed 
that if the ratio between the Holmberg radius of a galaxy 
and the half power beam width is smaller than ~ 7 to ~ 10, 
beam smearing can have serious effects on the steep inner 
parts of the curve. Rubin et al. (1989) make the point that 
if one wants to decompose the mass distribution of a galaxy 
into its disk and halo components one needs rotation curves 
of high accuracy in their inner portions, as "the maximum 
mass which can reside in the disk is constrained principally 
by the inner rise of the rotation curve." This statement is 
based on investigations into a sample of primarily steeply 
rising rotation curves, but remains true also for slowly ris- 
ing rotation curves, although, as we will show, to a much 
lesser extent. 

The rotation curves of LSB galaxies presented in 
BMH96 are not high-resolution rotation curves and beam 
smearing is a potential problem. The decrease in the slope 
of the rotation curves would disguise potentially steep rota- 
tion curves as gently rising solid-body curves. Here we will 
present some arguments why this effect is not significant in 
the BMH96 data, implying that the rotation curves of LSB 
galaxies are truly slowly rising. 



3.1 Evidence from BMH96 

The most direct evidence why beam smearing effects do not 
dominate the data from BMH96 comes from that data itself. 
If the data were severely affected by resolution effects, one 
would not expect to observe steeply rising rotation curves. 

Figure |l] shows the major-axis position-velocity dia- 
grams of LSB galaxies F579-V1 and F568-1. Both galaxies 
are at similar distances (77 vs 85 Mpc, a difference of less 
than 10 per cent); both were observed with the same tele- 
scope with identical beam sizes; both have an inclination of 
26 degrees, a comparable Hi distribution and almost identi- 
cal scale lengths. The effects of beam smearing are effectively 
identical in both galaxies. 

It is clear that the rotation curve of F579-V1 rises much 
more steeply than that of F568-1. The fact that we do ob- 
serve a steeply rising rotation curve in F579-V1, shows that 
the approximately solid body rotation observed in F568-1 is 
not caused by beam- smearing. This makes it unlikely that 
the solid-body rotation observed in other LSB galaxies from 
BMH96 is caused by observational effects. It should be noted 
that the solid-body rotation is also seen in those LSB galax- 
ies observed at higher resolutions (e.g. F563-1 and F583-1) 
and is in general found in other high resolution observations 
of very late-type galaxies (see e.g. IC 2574 in Martimbeau, 
Carignan & Roy 1994). 



3.2 Modeling beam smearing 

Another way to quantify the effects of beam smearing is 
to construct model galaxies with known properties and 
"observe" these at increasingly lower resolutions. As beam 
smearing affects the full two-dimensional velocity field of 
a galaxy, it is necessary to construct complete model data 
cubes, smooth these to lower resolutions using a Gaussian 
beam and then construct velocity fields and derive rotation 
curves. 



3.2.1 Constructing models 

We constructed a number of model data cubes using the 
task GALMOD in the Groningen Image Processing System 
GIPSY. This program distributes "Hi clouds" in a specified 
data cube using an input radial Hi distribution and rotation 
curve as distribution functions. The galaxy can be given 
any inclination or position angle. This model can then be 
"observed" at any desired spatial and velocity resolution. 
For simplicity we adopted a uniform, constant density Hi 
distribution for all models. We could have used any Hi dis- 
tribution, but decided to opt for the simplest distribution to 
isolate the effects of beam smearing. 

The magnitude of the beam smearing effects also de- 
pends on inclination. At high inclinations the beam will 
cover a larger part of the velocity field off the major axis 
than at lower inclination, thus "diluting" the major axis ro- 
tation velocity. For a fixed beam size the effects of beam 
smearing will still depend on inclination, with the worst ef- 
fects occuring in the highest-inclination galaxies. 

In order to mimic the BMH96 observations, we chose 
an inclination of 40 degrees for our model galaxies, equal 
to the average inclination of the 12 "F"-galaxies in BMH96. 
To cover the full range of known rotation curve shapes we 
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Figure 1. Major-axis position velocity diagrams of LSB galaxies F568-1 and F579-V1. Both galaxies are at approximately the same 
distance, have identical inclinations and were observed with the same telescope with identical beam sizes. It is clear that the rotation 
velocity in the inner parts of F579-V1 increases more rapidly than in F568-1. If these observations suffered from heavy beam-smearing 
this difference could never have been observed. Contour levels are -2(t (dashed contours) and 2,4,6,.. cr contours (full contours). 



adopted three difTerent input rotation curves. These are 
shown in Fig. ^ Radii and velocities are all expressed as 
fractions of the maximum radius and velocity of the input 
models -Rmax and Knax- These were identical for all three 
models. 

The first model has a constant rotation velocity at all 
radii, and is the perfect example of a flat rotation curve. 
This is an extreme, but not unrealistic representation of the 
steeply rising rotation curves found in early-type HSB galax- 
ies (e.g. UGC 2885 or NGC 2841 in the compilation by Rhee 
1996). 

The second curve is a simplified version of that pre- 
sented in FaU & Efstathiou (1980): 

vir) - y^^^^{^.- 

The constant d was given a value of 0.2_Rinax- This model 
rotation curve (which we will refer to as the "FE curve") 
resembles that of a normal late-type galaxy (see e.g. the 
high resolution rotation curves of NGC 1560, NGC 3198 or 
NGC 2403 in Rhee 1996). 

The third curve is another extreme case: pure solid body 
rotation with a still increasing rotation velocity at -Rmax- 
This solid-body rotation is found in dwarf galaxies (see e.g. 
the well defined high resolution curve of IC 2574 in Martim- 
beau et al. 1994 or NGC 3109 from Rhee 1996). 

The full three-dimensional high-resolution model 
data cubes constructed using these curves, were each 
smoothed with Gaussian beams with FWHM sizes of 
(1/10, 1/8, 1/6, 1/4, 1/2) X _Rmax. This obviously resulted for 
each of the three cases in data cubes where the radius of the 
disk measured 10, 8, 6, 4 and 2 beams respectively. Velocity 
fields were constructed from these cubes using the GIPSY task 



MOMENTS. Analogous to the procedure outlined in BMH96, 
the rotation curves of each of the smoothed data cubes were 
determined from major-axis position- velocity diagrams. 

As noted before, all models were computed assuming 
model galaxy inclinations of 40 degrees. The origtial and 
smoothed de-projected curves are presented in Fig. 0. A first 
glance shows that indeed the effects of beam smearing are 
most severe for the flat rotation curve which has the steep- 
est inner "slope." The beam smearing effects are negligible 
in the solid body curve, where the only effect is an appar- 
ent flattening at the outermost radii. The FE curve retains 
its shape and slope down to a resolution of 4 beams. The 
input and output curves only start to differ signiflcantly at 
resolutions of 2 beams. 



3.2.2 Comparison with observations 

We now compare the smoothed models with some of the 
observations from BMH96. We will use the rotation curves 
of F568-3 and F568-V1. These galaxies both have inclina- 
tions of 40 degrees and are thus directly comparable with 
the model curves. Both curves have a resolution of slightly 
less than 4 beams and should thus be compared to the 4- 
beam models, as is done in Fig. |^. The observed curves were 
scaled using the -Rmax and Knax given in BMH96. 

It is clear that the observed curves are best matched 
by the FE curve, although a completely solid body curve 
cannot be ruled out at 4 beams resolution. The only way the 
observations can be reconciled with a rotation curves that 
are flat down to a small radius, is at a resolution of 2 beams, 
that is, a factor of 2 worse than the present observations. At 
a resolution of 4 beams a flat curve can be ruled out. 

The question that remains to be answered then, is how 
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Figure 2. Model rotation curves illustrating the effects of beam-smearing. The heavy lines show the input rotation curves: flat (left 
panel), late-type-like (center), and solid-body (right). The beam-smeared versions of these curves are shown as the dashed and dotted 
lines. These lower resolution curves were derived using the full 3-dimensional data cubes. The resolution expressed in beams is 10 (dotted), 
8 (short dash), 6 (long dash), 4 (short dash-dot) and 2 (long dash-dot). It is clear that the effects of beam-smearing are most severe in 
the flat case (left panel), while the effects of beam-smearing are only small in the inner parts of the late- type and solid-body curves. In 
the latter two cases the shape of the curves do not change dramatically down to a resolution of 2 beams. The small bump visible in the 
FE-curves at i? ~ 1 is caused by the limited velocity resolution used in the models, in order to mimic the BMH96 observations. 



likely it is that LSB galaxies have FE-like rotation curves. 
First, as is shown by the studies of Casertano & van Gorkom 

(1989) and Broeils (1992) of high resolution Hi rotation 
curves, only those galaxies with maximum rotation velocities 
of more than 180 km s~^ {Mb < —19) have steep inner rota- 
tion curves. The steep inner curves are usually interpreted 
as the dynamical signature of a strong bulge component. 
Galaxies that rotate slower have gradually rising rotation 
curves, resembling the FE model curve. The LSB galaxies 
from BMH96 have maximum rotation velocities between 50 
and 120 km s~^, luminosities Mb between —16 and —18, 
and have only the faintest trace of a central condensation. 
So, if LSB galaxies are late-type galaxies, as all of the avail- 
able evidence suggests, then it is not surprising that they 
have slowly rising rotation curves. 

This is supported by work from Coradi & CapacioUo 

(1990) who measured the inner slopes of optical rotation 
curves of 139 galaxies of a large range of Hubble types 
and luminosities. They found a strong correlation between 
the magnitude of this slope, Hubble type and luminosity: 
75 per cent of the galaxies in their sample with luminosi- 
ties 19 < Mb < — 18 have slowly rising, non-flat rotation 
curve^, while all of the galaxies fainter than —18 have 

* Note that the steep slope of the F579-V1 rotation curve is 



slowly rising rotation curves. The result holds for the re- 
lation with Hubble types: 80 per cent of the galaxies with 
Hubble type Sd and later have slowly rising rotation curves. 



3.3 Evidence from NGC1560 

As LSB galaxies are late-type galaxies, a third possibil- 
ity to investigate the effects of beam-smearing is by tak- 
ing a high-resolution observation of a late-type galaxy and 
smooth it to lower resolution. We have selected NGC 1560 
for this purpose. This galaxy is a LSB dwarf galaxy (/xo(-B) = 
23.2 mag arcsec"^, for other properties see Table 1), for 
which a high-resolution rotation curve is available (Broeils 
1992). This curve is shown in the left panel of Fig. ^. The 
beam size for this WSRT observation was 13" x 14", and 
the extent of the rotation curve is 38 beams. Combined with 
the small distance of only 3 Mpc, we can be that sure beam- 
smearing is of no importance in this observation. 

We have used one-dimensional smoothing with a Gaus- 
sian filter to bring the total rotation curve and those of the 



consistent with these results. F579-V1 is the most luminous 
galaxy in the "F" -sample from BMH96, and in fact falls in the 
19 < Mb < -18 category. 
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Figure 3. Comparison of the model rotation curves with observed 
curves. The three Hght hnes show, from top to bottom, the rota- 
tion curves of the fiat model, the FE model and the solid-body 
model at a resolution of 4 beams. Super imposed are the rota- 
tion curves of F568-V1 (squares) and F568-3 (circles). The latter 
curves are consistent with the FE-model, although the shapes 
suggest that a slightly steeper solid-body model is also possible. 
The flat model can be excluded. 



stars and gas to lower resolutions. As the effects of one- 
dimensional smoothing of the rotation curve are less severe 
than full two-dimensional smoothing of the velocity field, 
we have smoothed the curves to a resolution of 3 beams. We 
found by trial and error that this corresponds most closely 
to a full two-dimensional smoothing with a resolution of 4 
beams. 

The smoothed curve is shown in the right panel of 
Fig. ^ The other curves shown in the figure will be discussed 
in Section 5.1. A comparison of the high- and low-resolution 
curves shows that the only effects this beam smearing ex- 
ercise has had on the curves is a loss of detail. The shapes 
and amplitudes of the curves have remained unchanged. This 
again shows that slowly rising rotation curves can be studied 
at much lower resolutions than rapidly rising ones without 
a dramatic loss of information. We conclude that the rota- 
tion curves of the LSB galaxies are intrinsically slowly rising. 
Beam smearing effects do not dominate the BMH96 data. 



4 DISK HALO DECOMPOSITION 

We now return to the disk-halo decompositions and subse- 
quent mass modeling. Before giving the results of the de- 
compositions, we first describe the mass components used. 

Stellar disk For the stellar disk the B and _R-band pho- 
tometry presented in de Blok et al. (1995), and McGaugh & 
Bothun (1994) was used. The rotation curve of the disk was 
computed following Casertano (1983) and Begeman (1987). 
The disk was assumed to have a vertical sech-squared dis- 
tribution with a scale height zo = h/6 (van der Kruit & 



Searle 1981). The exact value of zq does not infiuence the 
amplitude and shape of the rotation curve to any signifi- 
cant degree. No bulge component was included as there is 
no evidence that bulges contribute significantly to the mass 
distribution in most LSB galaxies (see de Blok et al. 1995, 
McGaugh, Bothun & Schombert 1995). Prior to the decom- 
position the light distribution curves were resampled to the 
Hi resolution. 

Isothermal dark halo We assume a spherical pseudo- 
isothermal halo with a density profile 



(1) 



where po is the central density of the halo, and Rc the core 
radius of the halo. This density profile results in the rotation 
curve 



Vhalo(i?) = X UnGpoR, 



RcJ 



(2) 



The asymptotic velocity of the halo, Voo, is given by 

14o = ^/^7vGpoR%. (3) 



To characterize this halo only two out of the three parame- 
ters (po, Rc, Voo) are needed, as equation (3) determines the 
value of the third parameter. The mass of the halo integrated 
out to radius r is 



Md(r) 



G 



■ A-KpoRc 



Rc arctan 



(4) 



Because {M/L)i, is unknown we present disk-halo de- 
compositions using three different assumptions for {M/L)i,: 

(i) maximum disk (Sec. 4.2), 

(ii) minimum disk (Sec. 4.3), and 
(in) "Bottema disk" (Sec. 4.4). 

Each of the rotation curves was fitted using a program 
which determines the best-fitting combination of Rc, Voo, 
{M/L)i,, using a least squares fitting. Any combination of 
these three parameters could be fit, or kept fixed at some 
initial value. The program uses as input the rotation curve of 
the gas, the total measured rotation curve, and the rotation 
curve of the stellar disk computed from the observed light 
distribution (i.e. with an (M/L)* = 1) (see Begeman 1987). 



4.1 HSB and LSB samples 

We will compare the properties of our LSB galaxies with 
those of mostly HSB galaxies with well-defined Hi rota- 
tion curves. The latter have been taken from the compi- 
lation by Broeils (1992) and references therein. This "HSB 
sample" does not consist of just disks obeying Freeman's 
Law (i.e. galaxies with /io — 21.6B-mag arcsec"^) (Freeman 
1970): a large range of central surface brightnesses is found, 
from Freeman disks to galaxies with surface brightnesses 
approaching those of the LSB galaxies. The mean surface 
brightness of the HSB sample {{po) — 22.2) is still 1.5 mag- 
nitudes brighter than that of the LSB sample {{po) ~ 23.7). 
Comparing the two samples gives a good impression of the 
change in properties of spiral galaxies over a large range in 
surface brightness. Properties of the HSB sample are given 
in the bottom panel of Table 1. 
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Figure 4. High and low resolution rotation curves of NGC 1560. The left panel shows the original high resolution observation from 
Broeils (1992). The curves measures 38 beams, and is well-resolved. The right panel shows the same data but smoothed to a resolution 
of only 3 beams. It is striking that the overall shape and amplitude of the curves have remained virtually unchanged. This is discussed 
more fully in Section 3.3. The dashed and dotted lines show the maximum disk decompositions for both cases. The dotted lines represent 
the rotation curve of the gas component, the dashed lines the rotation curve of the maximum stellar disk. The dash-dotted lines are the 
rotation curves of the halos as derived from this decomposition. There is almost no difference between the low- and high-resolution halo 
parameters. This is discussed more fully in Section 5.1. The fitting parameters are given in Table 6. 



4.2 Maximum disk fits 

The maximum disk hypothesis attributes as much of the 
observed rotation velocity to the steUar disk as possible. It 
therefore always yields an upper limit to the (M/L)^ in a 
galaxy. For HSB galaxies the maximum disk hypothesis can 
explain almost all of the observed rotation velocity in the 
inner parts with just the stellar disk. For LSB galaxies the 
situation is quite different. Figure |^ shows the maximum 
disk decompositions of their rotation curves. The resulting 
disk and halo parameters are summarized in Table 2, to- 
gether with the masses of the various components, as mea- 
sured within the outermost radius i?max of the rotation curve 
(Table 1). For those galaxies where both R and B photome- 
try were available we have done separate decompositions for 
both bands, which are compared in Fig. ^. The core radii 
and asymptotic velocities correspond very well to each other 
in both sets, as do the maximum disk (A//L)* ratios when 
the B — R colors of the galaxies are taken into account. To 
compare those LSB galaxies for which only i?-band data is 
available with _B-band literature data, we will convert the 
measured maximum disk (M/L_r)* ratios to B-band values 
assuming an average B — R color of 0.88 (de Blok et al. 1995, 
cf. Fig. H). Our results do not depend on this conversion as 
the halo is the most dominant component and its parame- 
ters are therefore not very sensitive to the precise value of 
(M/L)*. The maximum disk decompositions for the HSB 
samples have been taken from Broeils (1992) and are for 
convenience summarized in Table 4. 



4.3 Minimum disk 

Many of the rotation curves of LSB galaxies are best fit- 
ted by a minimum disk, that is, under the assumption that 
{M / L)i, = 0. We have calculated the halo parameters under 
this assumption (Table 3). We will compare them with the 
maximum disk parameters in Section 5.4, but already note 
that there is little change in the structural parameters of 
the halos, due to the DM dominance. Although minimum 
disk is not a realistic assumption for HSB galaxies, we have 
for comparison also made minimum disk fits to the rotation 
curves of the HSB sample (Table 5). Not all galaxies could 
be fitted: especially the earlier types (with bulges) have ro- 
tation curves that rise too steeply to be fitted by a simple 
isothermal model with a finite core radius, thus clearly show- 
ing the need for a compact luminous component. 



4.4 Bottema disk 

Bottema (1995) has measured the stellar velocity disper- 
sions of a sample of HSB spiral galaxies and found that 
the dispersion is related to the maximum rotation veloc- 
ity. For a Freeman disk with a color B — V = 0.7 he de- 
rived (M/Lb)* — 1.8. The stellar disks of the galaxies in 
his sample supply approximately 63% of the total measured 
velocity at 2.2h (where the rotation curve of an exponen- 
tial disk peaks). These results were derived under the ex- 
plicit assumption of Freeman's Law (which holds for most 
of Bottema's galaxies). It is obvious that Bottema's orig- 
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Figure 5. Maximum disk rotation curve decompositions of the final sample of LSB galaxies. The dotted lines represent the rotation 
curves of the gas; the long dashed line those of the scaled stellar disk; the short dashed lines the rotation curves of the halo. The full 
line represents the total model rotation curve. Error bars are based on a combination of profile width in the position-velocity diagrams 
(BMH96) and the asymmetries between the rotation curves of both sides of the galaxies. Shown are the R band decompositions, except 
for F563-V2 where we show the B band data due to lack of R data, and except for F571-V2 where the minimum disk solution is shown, 
as no optical photometry is available for that galaxy. 



inal prescription cannot hold in LSB galaxies: 63% of the 
total rotation velocity is in many cases already more than 
the maximum disk solutions allow. In Appendix A we give a 
more general expression for (M/Ls)* of a galactic exponen- 
tial disk derived by Bottema (1997). The resulting (M/Lb)* 
only depends on color, as the effects of luminosity and sur- 
face brightness cancel in this model. In the absence of data 
on stellar velocity dispersions of low or intermediate sur- 



face brightness galaxies, this is the best we can currently do 
taking into account the actual gravitational potential of the 
disk. In the following we will refer to this prescription as 
the "generalized Bottema-disk." Of course, when applying 
this generalized scheme to HSB galaxies, the original 63% 
criterion is retrieved. 

The (M/L)* values given in column (20) of Table 3 
show that the range of values implied by the Bottema 



© 0000 RAS, MNRAS 000, 000-000 



10 W.J.G. de Blok and S.S. McGaugh 

Table 3. Disk and halo parameters LSB galaxies 



Maximum disk R-band Maximum disk S-band 



Name 


(M\ 
\ L )^, 




Voo 


PO 




Ma ( 


M) 




Voo 


PO 


Mi, 




W 






(.4) 


(.5) 






{°) 


(.9) 


(10) 




(12) 


(13) 


F563-1 


6.3 


3.6 


109.0 


17.2 


0.962 


3.529 


9.0 


4.7 


106.8 


9.6 


1.374 


3.063 


F563-V2 


* 


* 










2.7 


3.8 


115.2 


17.2 


0.658 


1.456 


F565-V2 


3.1 


4.6 


61.8 


3.4 


0.103 


0.310 


3.6 


5.6 


70.7 


2.9 


0.097 


0.331 


F568-1 


2.0 


2.6 


126.4 


43.7 


0.394 


4.163 


3.0 


2.6 


124.3 


40.8 


0.587 


3.993 


F568-3 


1.8 


5.4 


136.5 


11.8 


0.429 


4.190 


2.6 


6.6 


147.4 


9.1 


0.532 


4.305 


F568-V1 


2.0 


3.8 


135.9 


23.6 


0.207 


5.929 


2.6 


4.2 


137.2 


19.5 


0.332 


5.824 


F571-8 


0.3° 


■ 7.9 


249.9 


18.2 


0.097 


9.934 


* 






* 




* 


F571-V1 


3.0 


4.4 


78.3 


5.9 


0.253 


1.279 


3.7 


5.5 


86.5 


4.6 


0.228 


1.379 


F574-1 


0.9 


10.2 


172.1 


5.3 


0.216 


3.681 


* 


* 


* 


* 






F583-1 


1.5 


3.4 


97.8 


15.5 


0.045 


2.236 


* 


* 


* 






* 


F583-4 


0.9 


6.0 


104.2 


5.6 


0.065 


0.962 




* 




* 






U0128 


3.0 


6.4 


136.3 


8.3 


1.043 


14.32 


4.0 


6.7 


125.9 


6.6 


2.254 


12.120 


U1230 


0.7 


4.5 


102.7 


9.5 


0.187 


6.912 


1.0 


5.1 


100.8 


7.2 


0.457 


6.48 




4.5 


3.2 


91.8 


15.0 


0.675 


4.467 


* 




* 


* 




* 


U5750 


7 


4.5 


83.4 


6.4 


0.086 


2 526 


* 


* 


* 




* 


* 




3.0 


4.6 


190.7 


32.0 


0.444 


7.966 


* 


* 




* 




* 








Minimum disk 








Bottema disk, B-band 




Name 


M„ 




Voo 


PO 




( 
V 




Rc 


Voo 


PO 




Af* 


(^ A\ 


(15j 


(,16j 


(17) 


(18) 


(19) 




(20j 


(21) 


(22) 


(23) 


(24) 


(25) 


F563-1 


0.385 


1.5 


118.6 


111.2 


5.042 




1.7 


1.8 


113.4 


75.4 


4.507 


0.203 


F563-V2 


0.321 


1.9 


127.6 


80.9 


2.479 




1.3 


2.8 


118.9 


34.4 


1.858 


0.365 


F565-V2 


0.084 


2.8 


59.1 


8.4 


0.399 




1.3 


3.4 


59.1 


5.6 


0.355 


0.038 


F568-1 


0.557 


2.0 


130.0 


82.2 


4.731 




1.5 


2.6 


128.6 


45.4 


4.313 


0.246 


F568-3 


0.394 


3.4 


131.2 


28.2 


4.742 




1.6 


5.2 


135.7 


12.5 


4.219 


0.308 


F568-V1 


0.344 


3.1 


134.3 


35.0 


6.153 




1.4 


3.7 


135.5 


25.5 


5.972 


0.128 


F571-8 


0.202 


7.0 


235.3 


21.0 


9.725 




* 


* 


* 


* 




* 


F571-V1 


0.164 


2.5 


79.9 


18.6 


1.640 




1.4 


3.3 


78.2 


10.7 


1.448 


0.103 


F571-V2 


0.016 


0.8 


51.1 


75.5 


0.156 




* 




* 


* 




* 


F574-1 


0.485 


7.5 


150.6 


7.5 


3.704 




1.5 


12.6 


188.9 


4.1 


3.502 


0.284 


F583-1 


0.243 


2.8 


96.5 


21.8 


2.322 




1.5 


3.4 


97.9 


15.3 


2.235 


0.036 


F583-4 


0.077 


4.7 


96.4 


7.8 


1.008 




1.5 


7.4 


112.3 


4.3 


0.905 


0.085 


U0128 


0.882 


4.0 


137.7 


21.7 


16.020 




1.5 


5.2 


135.7 


12.7 


14.90 


0.444 


U1230 


0.812 


3.7 


103.4 


14.2 


7.275 




1.1 


5.8 


104.1 


6.0 


6.697 


0.294 


U5005 


0.406 


2.2 


102.9 


40.1 


6.034 




1.5 


2.4 


97.2 


31.4 


5.337 


0.182 


U5750 


0.140 


4.0 


87.6 


9.1 


2.900 




1.5 


5.5 


84.2 


4.3 


2.386 


0.146 


U5999 


0.252 


4.1 


192.1 


40.8 


8.520 




1.5 


4.3 


191.0 


36.0 


8.209 


0.175 


Note: (M/L), 


t is in 


units of {Mq /Lq 


) in the respect 


Ive bands; Rc 


is in 


units of kpc, Vx 


, in units of km 


s~\ Po 


units of 10""^ 


M0pc 


-3; all 


masses are 


in units of 10^" 


Mq. 















a: (Af/L)buigo = 



disk is very small compared to the maximum disk values. 
This is due to the relatively small range in B — V colors 
shown by the galaxies. We have used Bottema's original 
"63% recipe" (as presented in Rhee 1996) for galaxies with 
Hq < 23 mag arcsec"^, while the generalized Bottema disk 
was used for galaxies with fj,Q fainter than this value, where 
63% of the total rotation velocity is already more than max- 
imum disk allows. In practice the generalized disk was thus 
used for the LSB galaxies, NGC 247 and the three DDO 
galaxies. For NGC 55, NGC 3109 and NGC 5585 the Bot- 
tema recipes resulted in (M/Lb)* values larger than those 
implied by maximum disk and we have assumed the max- 
imum disk values. In at least one case (NGC 3109) it can 
be shown that the discrepancy arises because of a large dif- 
ference between the optical (Fabry-Perot) and Hi rotation 



curves. If the Fabry-Perot curve is used the maximum disk 
{M/ Lb)* value becomes larger than the Bottema value. 



5 RESULTS OF DECOMPOSITIONS 

5.1 The influence of resolution on disk- halo 
decompositions 

As we have established in Section 3 that beam-smearing 
is not determining the shapes of the rotation curves, we 
should not expect the results of disk-halo decompositions of 
the LSB galaxies to be influenced significantly by resolution 
effects. We test this by again looking at the case of NGC 
1560 (Figure @). 

Even the high-resolution decompositions of this galaxy 
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Table 4. Maximum disk decompositions of HSB sample 
(Brocils 1992) 



Name 




Rc 




PO 


Md 






DD0154 


1.2 


2 


59 


15.3 


0.39 


0.01 


0.04 


DD0168 


1.1 


2.7 


98 


24.2 


0.21 


0.02 


0.03 


DDO170 


2.7 


2.3 


75 


19.2 


0.85 


0.03 


0.06 


N55 


0.7 


7.9 


146 


6.4 


1.5 


0.29 


0.13 


N247 


4 


7.3 


136 


6.54 


1.33 


0.95 


0.1 


N300 


2 


6.32 


132 


8.3 


1.74 


0.42 


0.16 


N801 


4 


74.3 


302 


0.3 


19.14 


23.83 


2.86 


N1003 


1.1 


9.7 


133 


3.5 


7.83 


0.82 


0.82 


N1560 


3.5 


6.8 


133 


7.2 


0.94 


0.12 


0.1 


N2403 


1.8 


6.6 


154 


10.2 


6.26 


1.39 


0.43 


N2841 


5.1 


21.7 


308 


3.7 


116.1 


23.97 


2.09 


N2903 


2.9 


3.2 


166 


51 


12.52 


4.48 


0.34 


N2998 


2.1 


24.8 


242 


1.8 


26.95 


13.79 


3 


N3109 


0.5 


8.7 


141 


4.9 


0.76 


0.04 


0.06 


N3198 


3.5 


7.6 


156 


7.8 


11.25 


3.11 


0.69 


N5033 


5.3 


5.9 


170 


15.2 


18.23 


9.84 


0.91 


N5533 


6.8 


34.6 


255 


1 


52.89 


30.05 


3.2 


N5585 


0.5 


1.8 


99 


56.9 


1.63 


0.06 


0.17 


N6503 


1.7 


2.5 


115 


38.6 


5.74 


0.83 


0.22 


N6674 


4 


119.5 


655 


0.6 


53.36 


20.92 


3.9 


N7331 


5.8 


103 


982 


1.7 


32.42 


12.95 


1.58 


U2259 


3.4 


6.1 


137 


9.5 


0.94 


0.35 


0.05 


U2885 


2.1 


44.9 


382 


1.3 


1.87 


37.37 


4.8 



Note: (M/La)^, is in units of (Mq/Lq g); Rc is in units of 
kpc, Voo in units of km s~^, po is in units of 10~^ Mqpc"^; 
all masses are in units of IO^^Mq. See Broeils (1992) and ref- 
erences therein. 




0.5 1 1.8 2 2.2 0.5 
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Figure 6. Comparison of the core radii Rq (left panel) and asymptotic velocities (middle panel) Voo of the halos, and the maximum 
disk mass-to-light ratio (M/L)* of the stellar disk (right panel) resulting from separate B and _R band decomposition of the rotation 
curves of the subset of LSB galaxies for which both R and B band photometry were available. The B and R band decompositions are 
consistent with each other. The full drawn lines are lines of equality. The dotted line in the right panel is the relation between {M / Lb)* 
and [M/Lu)* for B - R = 0.88. 



are not entirely unambiguous, showing that even in a max- 
imum disk decomposition it is difficult to entirely constrain 
the [M /L)i, value. The quoted maximum disk [M/ Lb)* val- 
ues range from 3.5 (Broeils 1992, Chapter 10) and 4.1 (Bege- 
man, Broeils & Sanders 1991) to 4.9 (Broeils 1992, Chapter 
5). The latter case is however an attempt to fit the data 
without a dark halo, and we will discard that value. To be 
consistent with Tables 1 and 4 we adopt (M/Lb)* = 3.5. 



This maximum disk decomposition is shown in the left panel 
of Fig. ^ Corresponding halo parameters, taken from Table 
3 are repeated in Table 6. 

As described in Sect. 3.3, the rotation curves of NGC 
1560 were smoothed from a resolution of 38 beams to 3 
beams. We have attempted a maximum disk decomposition 
on this low resolution data, and find that, although again 
there is some freedom in [M/L)*, the best-fitting value for 
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Table 5. Bottema disk and minimum disk decompositions HSB sample 





Bottema disk" 










Minimum disk*" 




Name 








PO 




A/* 


Rc 






Md 


DD0154 


1 1 




58.3 


16.3 


0.39 


0.01 


1 1 
1. ( 
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1 no 
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i.O 


O 
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0.04 
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24.28 
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2.0 
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110.8 
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14.34 
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11.8 


38.55 
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281.2 


6.61 
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1.0 
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Note: (M/Lg)i, is in units of (Mq/Lq s) in the respective bands; Rq is in units of kpc, Vc 

in units of km , po is in units of 10"'^ Mqpc~^; all masses are in units of IO^^Mq. 

a: The Bottema disk decompositions were taken from Rhee (1996). 

b: A dash indicates that a minimum disk fit was not possible. 

c: (M/Lb)* not determined because of unknown contribution of bulge. 

d: {M/Lg)-!, is larger than maximum disk (Af/L^)*. 



Table 6. Halo parameters of NGC 1560 





high res. 


low res. 


{M/Lb)* 


3.5 


4 




133 


121 


Rc 


6.8 


6.2 


po 


7.2 


5.9 



the low-resolution data is ~ 4. We can exclude {M/Lb)* = 
3 or 5 as possible solutions for this particular the 
corresponding low-resolution models systematically over- or 
under-estimate the rotation velocities at certain radii. The 
low-resolution halo parameters are given Table 6, and are 
comparable with the high-resolution values. 

The value {M /Lb)* = 4 for the low-resolution case was 
determined by simply scaling the stellar rotation curve to 
match the innermost point. This procedure is thus for the 
low resolution case equivalent to a formal maximum disk fit 
in the high-resolution case. Studying high-resolution curves 
in detail shows that this scaling to match the innermost 
point is usually what happens in practice when doing a for- 
mal maximum disk fit. 



5.2 Maximum disk mass-to-light ratios 

As is clear from Fig. the rotation curves of the stellar 
disks in the LSB galaxies (even in the maximum disk case) 
do not even approach the maximum rotation velocities of 
the observed rotation curves. This is in sharp contrast with 
HSB galaxies where the rotation curve of the stellar disk is 
usually able to explain most of the rotation velocity in the 
inner parts. 

This is also illustrated in the top panel of Fig. ^. This 
shows the ratio of the peak velocity of the maximum disk 
at 2.2h and the total observed rotation velocity at that ra- 
dius. At central surface brightnesses brighter than 23 B-mag 
arcsec"^, the maximum disk can account for 85% ± 10% of 
the total rotation velocity, but this rapidly drops to only 
40% or less for LSB galaxies. 

This increasing inability of the maximum disk to ex- 
plain the observed rotation velocities is illustrated in the 
middle panel of Fig. which shows the ratio of the peak 
velocity of the maximum stellar disk and the maximum ob- 
served rotation velocity. This ratio systematically changes 
from ~ 0.9 at high surface brightnesses (i.e. the disk defines 
the maximum rotation velocity and the halo "conspires" to 
keep the rotation curve flat) to values of ~ 0.3 for the LSB 
galaxies (i.e. the halo completely determines the dynamics 
of the galaxies). The bottom panel of Fig. ^ again shows the 
ratio of of the peak velocity of the maximum disk at 2.2h 
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Figure 7. Top panel: the ratio of tlie peak rotation velocity of the 
maximum disk and the observed total rotation velocity at 2.2h. 
This ratio is approximately constant at ~ 0.85 for galaxies with 
IJ,(i{B) < 23 mag arcsec"'^ and drops steeply to values of ~ 0.5 
for the LSB galaxies. Even at 2.2h the importance of the disk de- 
creases toward lower surface brightnesses, despite higher values of 
(M/L)i,. Middle panel: the ratio of the peak rotation velocity of 
the maximum disk and the observed total rotation velocity at the 
outermost measured point. This ratio also changes systematically 
with surface brightness, from ~ 1 for the highest surface bright- 
ness galaxies, where the disk determines the maximum rotation 
velocity, to ~ 0.3 for LSB galaxies, where the halo completely de- 
termines the dynamics of these galaxies. Bottom panel: the ratio 
of maximum disk peak velocity and observed velocity at 2.2h. At 
Vmax < 150 km s~^ a range of ratios is observed, from ~ 0.85 for 
HSB galaxies to ~ 0.4 for LSB galaxies. The lack of points with 
low ratios at Vmax > 150 km s~^ might be caused by selection 
effects: LSB galaxies with such large rotation velocities have not 
yet been mapped in Hi. 



and the total observed rotation velocity at that radius, but 
as a function of Vi^ax- 

The maximum disk (M/L)* values are shown in Fig. ^ 
as a function of fj,Q , Knax and B — V. The relation between 
(M/Lb)* and V^ax (top-left) shows a clear segregation be- 
tween galaxies of different surface brightnesses: at each value 
of Vlnax the HSB galaxies have the lowest value of (M/Lb)* 
while the LSB galaxies always show high {M/Lb)* values. 
Also shown is the relation {M/Lb)* oc L^^^ which was pre- 
viously derived for HSB galaxies by e.g. Ashman, Salucci 
& Persic (1993) and Broeils (1992). Clearly the conclu- 
sion that late-type (i.e. lower surface brightness) galaxies 
have smaller maximum disk stellar mass-to- light ratios than 



early-type (i.e. higher surface brightness) galaxies was, in 
the light of this new data, a selection effect. The maximum 
disk [M/Lb)* is not a simple function of L (or Knax) alone, 
but also depends on surface brightness. 

The absence of a relation between {M/Lb)* and fio 
(top-center) is a result of comparing galaxies with different 
global properties: at each Vmax there is a large range in fj,o 
(and consequently a large range in {M/Lb)*)- Comparing 
galaxies at fixed means comparing galaxies of different 
masses and luminosities. 

The surface brightness segregation effect at fixed Vmax 
cannot be an artificial effect caused by the lower resolution 
of the LSB observations. We have shown in Section 3 that 
beam smearing effects are potentially most severe in HSB 
galaxies. If beam-smearing were to have affected the LSB 
curves this would imply even higher values for {M/ L)*, thus 
makingthe difference between HSB and LSB galaxies seen 
in Fig. y even more dramatic (see also discussion in de Blok 
& McCaugh 1996). 

5.2.1 The maximum exponential disk 

The segregation in {M/L)* can be explained from the prop- 
erties of the exponential disk and the systematics of max- 
imum disk fitting. As LSB galaxies obey the same TuUy- 
Fisher relation as HSB galaxies (Zwaan et al. 1995), this 
requires LSB galaxies to have larger scale lengths at fixed 
Vmax, in order to still reach the required value of L. The 
maximum rotation velocity of an exponential disk is given 
by w (X \/aoh, where ao is the central mass surface density. 
In practice the value of ctq is not known, and the distribu- 
tion of the light is used to compute the rotation curve of 
the disk. This rotation curve of the fight has a maximum 
rotation velocity of 

vi oc Y^Eo/i, (5) 

where Eo is the central surface brightness in linear units. 
(In the following we will denote the rotation velocity of the 
disk with v, while the observed total rotation velocity wiU be 
indicated with V.) In the rotation curve fitting process the 
mass-to-light ratio {M/L)* is introduced to convert surface 
brightness to mass surface density and compute the disk 
rotation curve. Under the maximum disk hypothesis {M/L)* 
is scaled to its maximum possible value. The peak rotation 
velocity of the maximum disk is then given by 

Wm OC y/{M/L)*T.oh. 
From Eqs. (5) and (6) it is clear that 
/MX 

\L )* wi 

We will now compare a HSB and a LSB galaxy at iden- 
tical positions on the TuUy-Fisher relation. We find for the 
peak velocity of the rotation curve of the light in the HSB, 
vf, that {vf)^ oc {T,^h^), and for the peak velocity in the 
LSB, v}:, that {v]')'^ oc (E^/i"^), where Eq'" and h^'^ are the 
central surface brightnesses and scale lengths of the LSB and 
the HSB, respectively. As their luminosities L oc Eq/i^ are 

identical we find /i oc Eq ^ and we can write as 



(6) 



(7) 



T,oh^ oc 



^0 



H^O 



EL 

^0 



^0 
yH 
^0 



(«?)'• (8) 
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Figure 8. Top row: The maximum dislc _B-band stellar mass- to-light ratio plotted versus maximum rotation velocity (left panel), 
central surface brightness (middle panel), and average B — V color (right panel). The open circles denote galaxies with /io(-B) < 21.9 
mag arcsec"^, the triangles represent 21.9 < (Uo(-B) < 23.2 mag arcsec"^ galaxies and the asterisks galaxies with /Uo(-B) > 23.2 mag 
arcsec"^. In the left panel the dotted line shows the relation (M/Lg)* oc L^/^. The relation has been arbitrarily shifted to follow the 
filled HSB data points. In the right panel the steep upper dotted curve shows the colour-mass-to-light ratio track from Larson & Tinsley 
(1978) while the shallow lower dotted curve shows the tra<;k as derived by Jablonka & Arimoto (1992). Both curves have been shifted 
to go through log(M/Ls)* = 0.55 a,t B — V = 0.55 (see bottom-right panel). Bottom row: Idem, but now showing the B-band stellar 
mass-to-light ratio as derived assuming a Bottema disk. 



Using Eqs. (7) and (8) we express {M/L)\ in term of 
(M/L)?. 




For galaxies at identical positions on the Tully-Fisher 
relation, the ratio of their maximum disk (M/L)* values 
thus depends on the ratio of the surface brightnesses and 
the ratio of their peak rotation velocities. This latter ratio 
depends on the shape of the rotation curves or, more specif- 
ically, the change of shape with surface brightness at fixed 



Vinax. We have already seen that the rotation curves of LSB 
galaxies rise more slowly than those of HSB galaxies, so we 
expect the peak velocities of the niaxinmm disk of the LSB 
galaxies to be smaller than those of the HSB galaxies. 

In order to get (M/L)^' > {M/L)f, Eq. (9) requires 

that 




It is easy to show that if both LSB and HSB galaxies had 
flat rotation curves, we find that (M/L)^ > (M/L)^ always. 
In practice, at each Knax a large range of v{2.2h) /V (2.2h) 
is found (Fig. 7). For the current sample, where a typical 
HSB galaxy has /io ~ 21.5, while a typical LSB galaxy has 
^0 — 24, we find that condition (10) requires 11^/1;^ ^ 0.5 
on average. We find the largest range of v(2.2h) /V (2.2h) 
at Vmax — 100 km s~^ in the bottom panel of Fig. 7. 
Taking the ratio of the average value for the LSB galax- 



© 0000 RAS, MNRAS 000, 000-000 



Dark and visible matter in LSB galaxies 15 



ies v{2.2h)/V{2.2h) ~ 0.45 and that of the HSB galaxies 
v{2.2h)/V(2.2h) ~ 0.85, yields a value of vl;,/v^ ~ 0.5. So 
although there will be individual cases where condition (10) 
will not hold (as is shown by the overlap of HSB and LSB 
points in the top-left panel of Fig. 8), on average condi- 
tion (10) is always fulfilled, and we find that the intrinsic 
properties of the exponential disk combined with the max- 
imum disk recipe always yield larger values of (M/L)^ for 
LSB galaxies than for HSB galaxies. No assumptions have 
been made on the stellar populations, ages, or evolution of 
the galaxies. The larger fitted maximum disk stellar mass- 
to- light ratios for LSB galaxies at fixed Vinax, therefore do 
not need to reflect the true evolutionary stellar mass-to-light 
ratio of the disk. 



5.2.2 Population effects? 

But is it in fact at all possible to explain the systematically 
higher maximum disk (M/L)-^ ratios of the LSB galaxies 
at fixed Knax in a consistent way by invoking systematic 
changes in the stellar populations? The systematically larger 
ratios imply an increasingly more evolved population with 
decreasing surface brightness: population synthesis models 
(Larson & Tinsley 1978) show that (M/Lb)* of an evolv- 
ing stellar population is in general an increasing function of 
time, and hence of increasing B — V color. 

Two examples are shown in in the right panels of Fig. 
where the upper dotted line represents (M/Lb)* as a 
function of the B — V color of a stellar population with 
a dechning star formation rate (Larson & Tinsley 1978). 
The lower dotted line is the track as derived by Jablonka 
& Arimoto (1992) who have used a model which also takes 
into account chemical enrichment. 

To a certain degree the maximum disk results agree 
with the trend of increasing (M/Lb)* with redder B — V: 
the lower envelope of the distribution of the points increases 
with redder color. Uncertainties in the low-mass end of the 
IMF make it very hard to derive absolute [M/ Lb)* values 
from any population synthesis model. This uncertainty may 
therefore shift the model track up or down, while e.g. metal- 
licity effects may make the slope of the track more shallow 
(as the Jablonka & Arimoto track shows). However, the ob- 
served range at each color is too large to be explained by 
models with a constant IMF and where the star formation 
history is well-behaved and a smoothly varying function of 
Hubble type. 

The high maximum disk (M/L)* are furthermore in 
direct conflict with observational evidence on colors (Mc- 
Gaugh & Bothun 1994, de Blok et al. 1995), metallicities 
(McGaugh 1994) and gas-fractions (McGaugh & de Blok 
1997). As the stellar populations of LSB galaxies are in an 
earlier evolutionary state than HSB galaxies (van den Hoek 
et al. 1997) one would expect them to have lower values of 
{M/L)i,. The values of (M/Lb)* as derived under the maxi- 
mum disk hypothesis are therefore not representative for the 
stellar population in LSB galaxies. 

5.2.3 Disk dark matter? 

Trying to explain the maximum disk results in a physical 
way, while having a stellar population with a lower (Af/L)* 



than implied by the maximum disk values, requires an ad- 
ditional dark component in the disk which traces the stel- 
lar population. This component has to become increasingly 
important towards lower central surface brightnesses. This 
might be an invisible baryonic component like e.g. cold 
molecular gas. There has already been much speculation 
about this in the literature (e.g. Pfenniger, Combes & Mar- 
tinet 1994). LSB galaxies would be the ideal test cases for 
such theories as the major part of their disks should then 
be made up of this baryonic DM. Again, this component 
should be of increasing importance towards lower surface 
brightnesses. 

This creates the unattractive picture of two dominant 
but distinct and differently distributed dark components: 
one that has to provide the extra mass in the disk to make 
it satisfy the maximum disk boundary condition, and one 
that has to provide the surplus velocity in the outer parts 
to produce the observed rotation curves. These two dom- 
inant components then have to "conspire" to produce the 
observed rotation curve. Either way, DM remains of extreme 
importance in LSB galaxies. 

5.2.4 Non-maximum disk 

Another possibility is that disks are simply non-maximal. 
Part of the halo DM is wrongly attributed to the disk sim- 
ply because the disk does not dominate the dynamics in the 
inner parts. This could in principle affect the derived values 
for the structural parameters of the halo, but this is only 
of minor importance in LSB galaxies. They are DM domi- 
nated, irrespective of the value of (M/L)^, so that changing 
the mass of the disk does not really matter. It is, however, of 
crucial importance in HSB galaxies where the disk is dom- 
inant, and the halo parameters therefore sensitive to the 
precise value of (M/L)*. 

The maximum disk hypothesis, in combination with the 
extendedness of LSB galaxies thus produces a balancing act 
where LSB galaxies have to have large (A//L)* ratios, or 
extra dark components, which are contrary to what we can 
derive from other evidence (stellar populations, metallicity, 
gas fraction). Maximum disk is therefore not the preferred 
way of making mass models in LSB galaxies. 

5.3 Maximum disk halo parameters 

For the LSB galaxies the asymptotic velocity Voo is well cor- 
related with the observed maximum rotation velocity Vmax, 
showing that although many of the LSB curves are still rising 
in the outermost point, they are close to their true maximum 
velocity. We will use this observational parameter instead of 
the fltted Voo . Figure ^ summarizes the halo parameters for 
the HSB and LSB galaxies under the assumption of maxi- 
mum disk. The core radius Rc is rather constant, both as 
a function of Vmax and /xq. The group of points at large Rc 
belong to HSB galaxies that are very much dominated by 
luminous matter, and are maximum disk in the true sense of 
the word: their inner rotation curve can be explained com- 
pletely by that of the disk. To avoid having a hollow halo in 
these galaxies, the maximum disk (Af/L)* is usually slightly, 
but somewhat arbitrarily, lowered. The halo parameters are 
extremely sensitive to this. In most cases these galaxies have 
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Figure 9. Isothermal halo fitting parameters for maximum disk fits (left panels) and Bottema disk fits (right panels). The open circles 
represent the HSB sample, the asterisks the LSB sample, po is expressed in units of 10~^Mq pc~^; Rc in kpc; (tq in 10~^ Mq pc~-^; 
Vmax in km s~^; and /io(B) in mag arcsec"^. 



a bulge, that also has to be accommodated in the fitting pro- 
cedure, with its own uncertainties in its (M/L)*. The HSB 
points at large Rc are thus uncertain. 

The central halo densities po as a function of /xo tend to 
be somewhat higher than those of HSB galaxies, but again, 
the uncertainty in the decompositions of the HSB galaxies 
makes it difficult to say anything more definite. 

The ratio between core radius and optical scale length 
changes with surface brightness from Rc — Sh for the HSB 
galaxies to Rc — h for the LSB galaxies. In the maximum 
disk picture the optical disks of LSB galaxies are massive 
and extend further out into the halo. 



5.4 Minimum disk 

We can illustrate that the halo parameters derived for the 
LSB galaxies are robust values, by comparing the values de- 
rived for maximum disk and minimum disk. This is done in 
Fig. ^ where the core radii and central densities as derived 
using these two extreme hypotheses are compared. The dif- 
ference in maximum and minimum disk halo parameters is 
clearly a strong function of surface brightness. The largest 
difference is observed for the HSB galaxies, where the cen- 
tral density po changes by more than an order of magnitude. 
A similar conclusions applies to the core radius. This is in 
sharp contrast with the core radii derived for the LSB galax- 
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Figure 10. The central densities (top) and core radii (bottom) 
of the halos HSB and LSB galaxies under the maximum disk 
assumption (open circles) and minimum disk assumption (filled 
circles). The core radii and densities of the HSB galaxies are more 
affected by the respective assumptions than those of the LSB 
galaxies. Note the different slopes of the trends under the different 
hypotheses. The drawn lines are least-squares fits. 



ies. The minimum disk values differ from the maximum disk 
values by less than a factor of two. The slope of the trend 
with surface brightness changes from negative to positive 
when going from maximum disk to minimum disk. This is 
entirely caused by the large change in parameters of the HSB 
galaxies. 

The DM dominance makes the parameters of the LSB 
galaxy halos insensitive to the precise (M/L)* value. The 
values we derive for the halo structural parameters, but also 
for the dark mass in these halos, are therefore likely to be 
very close to their true values. 



5.5 Bottema disk mass-to-light ratios 

The most important property which distinguishes the Bot- 
tema disk from the maximum disk is its small range of 
(M/Lb)*- This is immediately apparent in Fig. ^. The Bot- 



tema disk typically implies values of (M/Lb)* between 1 
and 2. In general the reddest galaxies have the highest mass- 
to-light ratios. 

The striking systematic offset in [M / L)i, at fixed Vmax 
between HSB and LSB galaxies for the maximum disk hy- 
pothesis has disappeared (compare with top left panel in 
Fig. ^ . At fixed Knax galaxies now have approximately iden- 
tical (Af/Z/)* ratios. The majority of the data points has a 
much smaller spread and is in much closer agreement with 
stellar population models (Larson & Tinsley 1978; van den 
Hoek et al. 1997) than the maximum disk values. The need 
for an extra dark component in the disk has disappeared. 

5.6 Bottema disk halo parameters 

Whereas in the maximum disk case the galaxies with the 
highest values of Vmax have the lowest central halo densities, 
in the Bottema disk case they have the densest and most 
compact halos. The central densities of the galaxies with low 
Vmax are, as expected, less affected. The large change in halo 
parameters of HSB galaxies compared to their maximum 
disk values, results in trends of halo density and core radius 
opposite to those derived for maximum disk. 

Paradoxically, the uncertainties in the decompositions 
of the HSB galaxies makes it very hard to make definite 
statements on the halo properties and their possible changes 
with surface brightness. 

5.7 Halo parameters and galaxy evolution 

Generalizing (and running the risk of over-interpreting) the 
maximum disk and Bottema disk decompositions leads to 
two distinct different sets of relations between the dark and 
visible components of galaxies. The maximum disk decom- 
positions suggest that at fixed Vmax all halos are equal, while 
the Bottema disk decompositions suggest that the sizes and 
densities of DM halos are directly related to those of the 
optical disks that inhabit these halos. In summary: 

(i) Maximum disk Towards lower surface brightnesses 
halos become denser and smaller. Disks of lower surface 
brightness have higher (Af/L)* at fixed Vmax, and extend 
further out into the halo. 

(ii) Bottema disk Towards lower surface brightnesses 
halos become more diffuse and more extended. Disks of lower 
surface brightness have similar (Af/L)* at fixed Vmax and are 
found in more massive halos. 

Option (i) implies that the dimmest and most extended 
(LSB) galaxies live in the most compact halos. Presumably 
the collapse of their massive, dark disks was not as complete 
as in HSB galaxies (because of the larger angular momenta 
of LSB galaxies?). 

With option (ii) the properties of the baryonic part of 
the galaxy are reflected in those of the host halo: dim, diffuse 
and extended galaxies inhabit extended, low density halos. 
In both cosmologies the dominance of DM increases towards 
lower surface brightnesses. 



5.7.1 Maximum disk cosmology 

If, for the moment, we assume the maximum disk results 
to hold, what are the implications for the evolution of LSB 
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galaxies? In other words, what if the only difference between 
HSB and LSB galaxies at fixed Vmax is that LSB galaxies 
stick out further into otherwise identical halos? 

The picture which can then be sketched is deceivingly 
simple. For a HSB and a LSB galaxy at identical positions 
on the TF relation, the LSB galaxy has to be more ex- 
tended optically (Zwaan et al. 1995), and in order to have 
the same luminosity as the HSB galaxy, the amount of past 
and present star formation per area (the surface brightness) 

has to be lowerQ. The LSB optical disk simply encompasses 
more DM than the HSB optical disk, due to its larger extent 
in an otherwise identical halo, resulting in a LSB optically 
extended galaxy with a high total (i.e. DM and luminous 
matter) M/L. The low star formation rate (as averaged over 
the lifetime of the galaxy) implies slow evolution in the LSB 
galaxy and a low (M/L)* 0, which is difficult to reconcile 
with the maximum disk hypothesis assumed to hold to de- 
rive the these conclusions. 

Current ideas on the formation and evolution of LSB 
galaxies (Bothun et al. 1993, McGaugh 1992, Mo et al 1994, 
Dalcanton et al. 1997) suggest that LSB galaxies form from 
lower amplitude density peaks than HSB galaxies. As they 
are also found to be more isolated from their neighbors, and 
to have a lack of companions with respect to HSB galaxies, 
they have presumably also undergone less interaction with 
other galaxies. The extended gas disks also suggest that the 
collapse of the baryonic matter has been much less efficient 
than in actively star forming HSB galaxies. 

Although the Bottema disk solution may seems the 
most one, both solutions do have their problems. The maxi- 
mum disk solution has to explain why galaxies in otherwise 
identical halos have different collapse-histories as a function 
of surface brightness, while the other solution has to ex- 
plain how galaxies with different halo masses but identical 
luminosities conspire to still end up at the same position on 
Tully-Fisher relation. 

One conclusion is firm however, the dynamics of LSB 
galaxies are fundamentally different from those of HSB 
galaxies. For example, the rotation curves of LSB galaxies 
can not be explained in fio = 1 cosmological simulations 
(Moore 1994, Navarro 1997), which is a problem, certainly 
in view of the implied numerical richness of LSB galaxies 
(McGaugh 1995). 



6 MASS RATIOS 

6.1 Importance of atomic gas 

The increasing dynamical importance of the gas component 
can readily be seen in Fig [ll| a and b. Galaxies of lower sur- 
face brightnesses have larger values of Mgas/M*. A similar 
trend with Hubble type was already noticed by e.g. Broeils 

t The bluer colors of LSB galaxies are explained by metallicity 
effects and the less prominent old population. An equal amount 
of star formation has a much larger effect on the colors of a LSB 
galaxy than on a HSB galaxy (de Blok et al. 1995). 
■t Here we specifically do not include Hi and other potentially star 
forming material; the low (M/L)* simply reflects that the LSB 
has not had enough time to build a large old stellar population 
as in HSB galaxies. 



(1992). The trend is clear in both the maximum disk and 
Bottema cases, although of course the gas-to-stellar mass 
ratios derived using maximum disk stellar masses tend to 
be larger. The ratio ranges from ~ 0.2 for the HSB galaxies, 
to ~ 2 in the LSB galaxies. Clearly in LSB galaxies the gas 
component can be of a larger dynamical importance than 
the stellar component, even in the maximum disk case. This 
effect is more pronounced for the Bottema decompositions 
because of the limited range in (M/L)*. A fuller discussion 
of the gas fraction over a large range in surface brightness 
is presented in McGaugh & de Blok (1997). 

Mgas/Afdark (as computed within the outermost mea- 
sured point of the rotation curve, i?max (see Table 1), and 
assuming Bottema disk) is shown in Fig. pdj c and e. Pan- 
els d and f show Mgas/Mdark measured within 5/i. There is 
clear trend of increasing gas-richness towards low values of 
Knax, an effect which is not entirely cancelled when mea- 
sured within 5/i. Slower rotating galaxies contain relatively 
more gas. The large range in /io found at each Vmax ensures 
that no clear trend is visible of Mgas/Mdark with /^q. 

6.2 Baryon fraction and dark matter fraction 

The total luminous mass Afium is defined as the sum of stel- 
lar mass and the gas mass. We here take the stellar mass 
as implied by the Bottema disk. The ratio of the dark and 
luminous mass is related to the maximum velocity (or lumi- 
nosity) (Persic & Salluci 1991, Broeils 1992). Galaxies with 
low maximum velocities and low luminosities are more DM 
dominated than high luminosity galaxies. This effect is not 
very pronounced when Mdark is computed within the outer- 
most radius (Fig. and i): the extended and well-observed 
rotation curves of HSB galaxies encompass a relatively much 
larger fraction of the halo. The effect becomes much more 
pronounced when measured within a fixed number of scale 
lengths (Fig. |l^h and j). 

The extreme DM domination of the lowest surface 
brightness galaxies has more general implications. One ex- 
ample is the baryon fraction ft of the universe. This is of di- 
rect cosmological significance since Slo = ^b/ ft, where fit, is 
the baryon density of the universe: a measure of /t, combined 
with fii, from primordial nucleosynthesis directly yields SIq. 

A value of /(, thought to be representative of the univer- 
sal baryon fraction is given by X-ray observations of galaxy 
clusters (White et al. 1993; David et al. 1995; White & 
Fabian 1995): 0.1 < ft < 0.2. A similar result is obtained for 
compact groups of galaxies (Pildis et al. 1995). Primordial 
nucleosynthesis gives fib = 0.022 (e.g. Walker et al. 1991); 
together these suggest flo ~ 0.2. 

Now consider the rotation curve data. Dark matter is 
required to keep the rotation curves flat, and since there is 
no indication of declining rotation curves, the halos must be 
substantially larger than the last measured point. The dark 
to luminous mass within -Rmax therefore gives an upper limit 
on ft (Fig. 0), through ft < (Mium/Mtot) and therefore 

To be conservative, let us consider the maximum disk 
case. For most HSB galaxies, iZmax is not sufficiently large to 
provide an interesting limit. However, the DM domination 
of LSB galaxies improves this situation. In our sample, there 
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are a few galaxies for which the upper limit on fb is lower 
than the lower bound of the range of the cluster baryon frac- 
tions. These are F568-V1, UGC 5750, UGC 5999 (all with 
ft < 0.08), and F571-8 {ft < 0.05) (although this latter 
value is uncertain due to the inclination of the galaxy). The 
rotation curves of all four of these galaxies are flat or rising 
at the last measured point, so the dark halos must extend 
well beyond J?max to avoid having a perceptible effect. The 
baryon fractions of these galaxies really are different from 
those derived for clusters of galaxies. 

It is not obvious why the baryon fraction should differ 
from halo to halo. If for the moment we assume our limit to 
be representative of the universal fraction, then the limits 
on fib and fb give Qo > 0.45. However, the methods used 
to obtain the cluster results have now been thoroughly an- 
alyzed by many groups (e.g. White & Fabian 1995; Evrard 
et al. 1996) and seem robust. 

It may be that really varies from halo to halo, per- 
haps depending on scale. What physics might cause this is 
unclear. Dissipation acts in the wrong direction: if galaxies 



have dissipated a lot, their halos are very large, worsening 
the problem. It has been suggested that low mass galax- 
ies experience significant wind driven mass loss (e.g. Dekel 
& Silk 1986) which might rid them of some baryons. How- 
ever, this clearly has not occurred in the galaxies wc are 
discussing. Their gas has not been swept out; they are very 
gas rich and massive. They show no evidence of blow-out 
events (Bothun et al. 1994). It is therefore entirely ad- hoc 
to vary fb, and this would add yet another free parameter 
to the current cosmological models. 

Of course, until the nature of the dark matter is known, 
it is conceivable that some or all of it is baryonic. Thus, an- 
other possibility is that some fraction of halos is baryonic, 
and relatively more baryons have been incorporated into the 
halos of LSB galaxies, perhaps as MACHOs. How and why 
this should occur is again inobvious. LSB galaxies show lit- 
tle evidence of an early epoch of star formation which ought 
to provide the remnants that become MACHOs (McGaugh 
& Bothun 1994; de Blok et al. 1995), and it is not obvi- 
ous that these should be distributed in a halo. Completely 
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Figure 12. Top panel: Average volume density of LSB galaxies as 
measured within 5 disk scale lengths. Bottom panel: Central mass 
surface density implied by the maximum disk fits. Both quantities 
depend on surface brightness. 



non-baryonic matter is still required unless primordial nu- 
cleosynthesis or large scale estimates of Qq are wrong. The 
limits derived from the LSB galaxies are, however, hard lim- 
its, and they should be reconciled with the cluster values. If 
that does not prove to be possible, this might at last resort 
even indicate the break down of the dark matter hypothesis 
(Sanders 1997). 



7 SURFACE AND VOLUME DENSITIES 
T.l Total mass volume densities 

A question which remains to be answered is whether the low 
surface densities now observed in LSB galaxies also imply 
that they are truly low volume density objects. This question 
is extensively discussed in de Blok & McGaugh (1996) where 
the LSB galaxy UGC 128 and HSB galaxy NGC 2403 (which 
lie at identical positions on the TF relation) are compared. 
The conclusion reached there is that the total mass volume 
density enclosed within a fixed number of scale lengths is a 
factor ~ 10 lower in UGC 128 than in NGC 2403. We refer 
to that work for a discussion on the preferred length scale 
for measuring densities in galaxies. 

We have extended this analysis to our samples of HSB 
and LSB galaxies. The total mass density within 5h was 
computed and shown as a function of surface brightness in 
Fig |l^. The conclusion reached in de Blok & McGaugh is 
confirmed: the difference in (p) between Freeman like galax- 
ies and LSB galaxies with ~ 24 is again a factor of 10. 

Although a large range in mass density is inferred at 
each value of ^o, LSB galaxies are "diffuse" objects, even in 
the maximum disk case. The increase in optical scale length 
towards lower surface brightnesses is much faster than any 



change in the DM length scale. The increase in area spanned 
by the optical disk can therefore not be compensated by a 
corresponding increase in total matter enclosed by the op- 
tical disk. The enclosed mass increases with V^h, while the 
enclosed area (volume) increases with (h^)- The increas- 
ing scale length therefore ensures that for LSB galaxies the 
average volume or surface density enclosed by the disk is 
always lower. 



7.2 Maximum disk baryonic surface densities 

The maximum disk results are a good way to demonstrate 
that LSB galaxies have truly low matter surface density 
disks. That is, they are not normal surface density disks 
that happen to have less light per square parsec, but the 
underlying mass surface density distribution is more diffuse 
too (de Blok & McGaugh 1996). 

In Fig |l^ we show the central surface densities calcu- 
lated by multiplying the central surface brightness with the 
respective (M/L)* ratios. The resulting surface density is 
the maximum density the disk can have within the con- 
straints of the observed rotation curve. These values have 
not been corrected for the presence of Hi, as the surface 
density of the Hi gas becomes progressively lower towards 
lower surface brightnesses (BMH96). The increase in maxi- 
mum disk (M/I/)* ratio is not fast enough to compensate for 
the decrease in surface brightness: an effect which will only 
become more pronounced when hypotheses other than max- 
imum disk are used. Maximum disk is the most conservative 
scenario, but already here it is clear that the mass surface 
density in the disk is a function of surface brightness. 



8 SUMMARY 

We have presented rotation curve decompositions of a sam- 
ple of LSB galaxies and compared the properties of their 
constituent mass components with those of a sample of HSB 
galaxies. 

LSB galaxies are very much dominated by dark mat- 
ter thus proving to be the ideal test-cases for theories on 
the structure and formation of dark matter halos. The halo 
parameters derived for these galaxies are insensitive to the 
assumed mass-to-light ratio of the stellar disk, ensuring that 
the value we derive are in all probability close to the true 
values. 

This is in sharp contrast with HSB galaxies, where the 
halo parameters are sensitive to the assumed stellar mass-to- 
light ratio. Comparison of the HSB and LSB samples shows 
that the maximum disk solutions imply an increasing stel- 
lar mass-to-light ratio for the stellar disk and an increasing 
compactness for the halo as the central surface brightness 
decreases, while the Bottema disk solutions imply exactly 
the opposite. We argue, based on what is known about the 
evolutionary state of LSB galaxies from other sources, that 
the maximum disk solution is not representative. To explain 
the maximum disk results in a physical way, one has to as- 
sume large amounts of baryonic (dark) matter in the LSB 
disks. This would result in massive, evolved disks, which is 
not consistent with the observed past and present low star 
formation rates of LSB galaxies. The Bottema disk solutions 
give lower (M/L)* values, more consistent with what can be 
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derived from other evolutionary probes (metallicity, gas frac- 
tion) . They imply that LSB galaxies live in extended diffuse 
halos. The halos of LSB galaxies are thus fundamentally dif- 
ferent from those of HSB galaxies of the same mass. This is 
a challenge for cosmological N-body simulations which will 
have to address this problem, especially considering the nu- 
merical richness of LSB galaxies. 

LSB galaxies are thus slowly evolving, isolated, diffuse, 
low-density and extremely dark matter dominated galaxies. 
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APPENDIX A: THE GENERALIZED 
BOTTEMA-DISK 

Al Introduction 

Bottema (1995) has shown from measurements of the stel- 
lar velocity dispersions in a sample of HSB galaxies that the 
maximum rotation velocity of the stellar disk (i.e. the ve- 
locity occuring at 2.2h) is on average equal to 63% of the 
measured rotation velocity at that radius. 

This result was derived under the explicit assumption 
of Freeman's Law (which applies to most of the galaxies in 
Bottema's sample). It is therefore not possible to directly 
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apply Bottema's result to our sample of LSB galaxies, es- 
pecially if one realizes that 63% of the observed rotation 
velocity in many LSB galaxies exceeds the maximum disk 
solution. 

In Bottcma (1997) a rc-analysis and generalization of 
these results is presented, also taking into account surface 
brightness and color effects. We here summarize the main 
points of that analysis. 

A 1.1 Poor man's population synthesis 

The disk is assumed to be unaffected by dust and contains 
only two populations (van dcr Kruit 1986; Bottcma 1988): 
a young component with a color {B — Vy = —0.03 and an 
old component with {B — V)° = 0.97. The old component 
furthermore contains all of the disk mass. 

It is straightforward to show that this implies that Sy = 
0.9735'Ij and S'^ = 2.4451;, where are the fluxes of the 

old and young populations in the B and V bands. 

The absolute magnitude of the old disk in the B-band, 
Mg''', is related to the total absolute magnitude M^"' by 

M^"" = MT - 2.5 log(S'|/5'B*) (Al) 



It is straightforward to show that this equation and the 
colors of the populations imply that 

Mi'a = Mb- 2.5 log ( ^7 °//^* ) = - cf, (A2) 
where A* is defined as iq-o-^^b-v) ^ 



A 1.2 The exponential disk 

For an exponential, locally isothermal, constant M/L disk 
with density distribution p{R,z) = p{R,0) sech^{z/zo), and 



2 - / 

and 



(A3) 



(A4) 



where (u^) is the central velocity dispersion, and uo. So 
and Zo are the central surface density, surface brightness and 
scale height, respectively. 

For an infinitely thin stellar disk the maximum rotation 
velocity is given by 



Wmax = 0.881/ nGaoh. 

Combining Eqs. (A4) and (A5) yields 



(A5) 



(A6) 



As the old disk contains all the mass Vmax.od = ^^max- 
(M/L)od is assumed to be the the same for all old disks. 

A2 Bottema disk for Freeman galaxies 

For Freeman galaxies /io is constant. The color is also as- 
sumed to be constant aX B — V = 0.7. Eq. (A7) then yields 
Mod = -10 log(t;inax) + C, or 

t,^3^ = 10°i^l0-°i^od. (AS) 
Combining Eqs. (A6) and (AS) yields 



{vi) 



2\5 



R=0 



= A- 



10 



-O.lMod 



where A is defined as 
^ = O.SS- lO""-^"^. 



(A9) 



(AlO) 



Following Bottema (1995) the empirical relation between lu- 
minosity and scale length ratios is adopted h/zo — 0.6 Mod + 
17.5 together with the corresponding values A = 0.75, 
C = 0.69 and (M/Lb) = 1.85. This results in 



-O.lMod 



(All) 



A3 General Bottema disk for non-Freeman 

galaxies 

The result derived above strictly speaking only applies to 
galaxies that obey Freeman's Law. To generalize to different 
surface brightnesses rewrite Equation (A7) as 



''max = O.SttG^So.f 



Sod 




if).,'-'"'' 



(A12) 



where Eo,f is the constant Freeman central surface bright- 
ness of 136 Lr.) pc^^. The scale lengths of the old and young 
populations are assumed to be equal, so that Sg'^/Eo = 
Lod/itot = 10°-*=-^. This results in 



= const 



10' 



,0.1c/ 



10 



-0.1M„d 



(A13) 



For Bottema's sample, with central surface brightnesses 
approximately equal to the Freeman value and colors close 
to B — V = 0.7, an average correction factor (c/) ~ —0.5 
can be defined. Combination of Eqs (A6) and (A13) then 
yields 



2.i _ const 



So,. 



4 lQO-l(cf+{cf)-M„^) 



(A14) 



Again using Bottema's relation h/zo = 0.6Mod -|- 17.5 
then implies that 

const X lO°-^<''f'> = 1.17 



A1.3 The old-disk Tully-Fisher relation 

Equation (A5) and the identity L = 2n'Eoh^, where L is the 
total luminosity, Eo the central surface brightness in linear 
units, and h the scale length, imply an old disk Tully-Fisher 

relation 



0.37rG 



1j / od 



(A7) 



For the maximum rotation velocity and velocity disper- 
sion of a galactic disk one then finds (using Eqs. (A13) and 
(A14)) 



Vmax = 1.17 

and 



-'0,F 



,0.1{c/-<c/)) j^Q-().lM„d 



(A15) 
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{vl)i^^ = 1.33,/5 ( ^\ ' ioO i(=/-<=«-^od) (A16) 
V h, \ijo,f J 

Conversion to observables using Eq (A2) then yields 



So 



3 



W = 1.17 1^1 ■ l0°-2(^^+°-2^) io-o i<- , (A17) 

A4 Mass-to-light ratio 

From Eqs. (A4) and (A16) one can derive that 

(^) = 28.1 . lO-°-2^«.0 . 10°-*=^ io-o-2<=/> , (A18) 
With Mb,q = 5.48 one finds 
^"l = 1.79-10°-*=^10°-^(°-«-<=^>' =2.84.10°-*=^ 

B 

= 1.936 • lO°-*(^-^' - 1.883. (A19) 



L 



The mass-to-light ratio does not depend on the surface 
brightness of the disk, which is due to our assumption that 
the old disk population has the same mass-to-light ratio for 
all disks. Disks with identical colors therefore have identical 
mciss-to-light ratios, irrespective of surface brightness. 
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